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ABSTRACT 

We study halo clustering bias with second- and third-order statistics of halo and mat¬ 
ter density fields in the MICE Grand Challenge simulation. We verify that two-point 
correlations deliver reliable estimates of the linear bias parameters at large scales, 
while estimations from the variance can be significantly affected by non-linear and 
possibly non-local contributions to the bias function. Combining three-point auto- 
and cross-correlations we find, for the first time in configuration space, evidence for 
the presence of such non-local contributions. These contributions are consistent with 
predicted second-order non-local effects on the bias functions originating from the 
dark matter tidal field. Samples of massive haloes show indications of bias (local or 
non-local) beyond second order. Ignoring non-local bias causes 20 — 30% and 5 — 10% 
overestimation of the linear bias from three-point auto- and cross-correlations respect¬ 
ively. We study two third-order bias estimators which are not affected by second-order 
non-local contributions. One is a combination of three-point auto- and cross- correl¬ 
ation. The other is a combination of third-order one- and two-point cumulants. Both 
methods deliver accurate bias estimations of the linear bias. Furthermore their estima¬ 
tions of second-order bias agree mutually. Ignoring non-local bias causes higher values 
of the second-order bias from three-point correlations. Our results demonstrate that 
third-order statistics can be employed for breaking the growth-bias degeneracy. 

Key words: linear bias, non-linear bias, non-local bias, clustering, higher order auto- 
and cross-correlation functions 


1 INTRODUCTION 

With the increasing amount of data coming from current 
and future large-scale galaxy surveys, errors on the ob¬ 
served statistical properties of the spatial galaxy distribu¬ 
tion are rapidly decreasing. This high level of precision re¬ 
quires at least the same level of accuracy in the modelling 
of the corresponding observables. An important observable 
is the growth of large-scale density fluctuations with time, 
which is sensitive to the universal matter density, the ex¬ 
pansion of space as well as to the gravitational interac¬ 
tion of matter at large scales. Measurements of the this 
growth therefore provide constrains on cosmological para- 


meters (e.s. iRoss et al, 

I2OO7I: ICabre and Gaztanaeal l2009l: 

Song and Percivall 200£ 

: Samushia et al.l 20121: Reid et al.l 

2012l:lde la Torre et al.ll2013l). Dossible deviations from Gen- 

eral Relativity (Gaztanaga and Lobdl200ll: iLue et al.ll2004l) 


or on alternative phenomenological description for the 
accelerated expansion, such as the effect ive field theory 
dSteigerwald et al]l2014l : IPiazza et al.ll2014l ). Growth meas¬ 
urements can be undertaken by comparing the second-order 


correlations ^ of galaxy distributions at different redshifts. A 
critical aspect of this approach is the bias between the cor¬ 
relations of galaxies and those of the full matter density field. 


This bias can eit her be predicted with the peak-background 
split model (e.g. Bardeen et al.l[l9^ : ICole and Kaiserlll989l : 
ISheth and Tor menll 19991 ). or directly determined from obser¬ 
vations using weak lensing observables, redshift space distor¬ 
tions or reduced third-order correlations at large scales. The 
latter method relies on the fact that such third-order cor¬ 
relations are independent of the growth at large scales, but 
sensitive to the bias. Third-order galaxy correlations there¬ 
fore have the potential to tighten constrains on cosmological 
mode ls from observations (e.g. iMarf 3 I2OIII : iMarfn et al.l 
I2OI3I) . However, how useful third-order correlation are for 
this purpose depends on the accuracy a nd the precision 
of the bias estimations they deliver (e.g. IWu et aD l2010l : 
lEriksen and Gaztanagall2015l ). 


The present work is part of a series of papers with 
which we aim to characterise and understand differences 
between various bias estimations using the huge volume 
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of the MICE Grand challenge simulation. We thereby as¬ 
sume that galaxies can be associated with dark matter 
haloes. This analysis includes bias measurements derived 
from a direct comparison between halo and matter fluctu¬ 
ations (Bel, Hoffmann, Gaztanaga, in preparation) as well 
as bias predictions from the peak-background split model 
llHoffmann et al.ll2015al l. Continuing this series we now aim 
at deepening our previous work on t he bias from third-order 
correlations (iHoffmann et al.ll2015bll . 

In this previous work we investigated growth measure¬ 
ments based on bias estimations from third-order halo auto¬ 
correlations (halo-halo-halo). Despite the positive and en¬ 
couraging outcome of this analysis, we pointed out some 
discrepancies between the first- and second-order bias es¬ 
timations derived from two different third-order statist¬ 
ics. The first third-order statis tics i s the reduced three- 
point correlation function Q (iFrvI Il984bh. The second 
one is based on a par ticular combination (ISzapudil 1 19981 : 
iBel and Mnxinonj 20121') of the sk ewness S 3 and the correl- 
ator C 12 (iBernardeau et al.| l2003i (referred to as r method, 
iBel and Marinonil l2012f i. which correspond to the three- 
point correlation for triangles with one and two collapsed 
legs respectively. We found that linear bias estimations 
from the three-point auto correlations over-estimate the 
true linear bias (probed by two-point correlations at large 
scales) by 20-30%. The r method delivers accurate lin¬ 
ear bias estimations, but with lower precision. Further¬ 
more, its bias estimation become unreliable when samples 
of massive haloes (M > Mq) are considered. Under¬ 

standing these discrepancies between different bias estimat¬ 
ors is crucial for constraining cosmological models with ob¬ 
served t hird-order galaxy statistics. In the light of previous 
studies (iM arieraandWzaatanasa l201l|; I Pollack et al.l I2OI2I : 
iGhan et al.ll2012l : iBaldauf et al.ll2012h we pointed out that 
such discrepancies might originate from non-linear and/or 
non-local effects on the different bias estimators. Hence, in 
the present paper, we focus on the analysis of third-order 
cross-correlations (halo-matter-matter) to decrease the im¬ 
pact of non-linearities on our bias estimators in order to in¬ 
vestigate the extension of the local bias model to a possible 
non-local component. 

We recall th e bia s expansion, introduced by 
iFrv and Gaztanagal lll993h . which assumes a local re¬ 
lation between the density contrast of matter and haloes 


5,, = F[5m] 


( 1 ) 


In the above relation we assume that 5h and Sm 
are coarse grain quantit ies obtained by appl yi ng a 
smoothing process (see iManera and Gaztah^al I2OIII : 
iBuchalter and Kamionkowskil ~999l . for a discusion of the 
effect of the smoothing). Inaccuracies of this deterministic 
relation (equation[T]) might arise from tidal forces in the mat¬ 
ter held, leading to a non-local contribution in the biasing 
relation. At second order it can be expressed as 


5,(x) = 61 |5^(x) + - (Sf)} + ■^e2 (x)|, (2) 

where 72 repres e nts the non - local bias parameter 
llGhan et al. I I2OI2I: iBaldauf et all l2012lb This non-local 


component depends on the divergence dv of the normalised 
velocity held {v/T-Llf) 

02(x) = - J /3i20^(qi)6'«(q2)IU[gi2J?]e*'*^^ "‘d^qid®q2, 

( 3 ) 

where /3i2 = 1 — ^ ) represents the mode-coupling 

between density oscillations with wave vectors q, which de¬ 
scribe tidal forces. W[qi 2 li\ is the Fourier transform of a 
spherical Top-hat window with radius R. In order to be con¬ 
sistent with the dehnition of second-order bias parameter, 
in this paper we shall refer to the non-local component of 
the biasing relation ([2)) using the quantity 32 = 2^. 

Evidence for signihcant contributions of such a non¬ 
local component to bias function has been reported in 
Fourier space for d ifferent simulations dGhan et al.l I2OI2I : 
IBaldauf et al.l I2OI2I I. However, it remains unclear how 
strongly these non-local contributions affect the bias and 
consequently third-order statistics of large-scale halo distri¬ 
butions in conhguration space. We address this latter ques¬ 
tion in the present study and suggest possibilities to employ 
third-order statistics for accurate bias measurements, inde¬ 
pendently of non-local bias. 

This paper is organised as follows. In Section [2] we 
briefly present the simulations on which our analysis is based 
on, in Section [3] we present the bias estimators implemen¬ 
ted in this paper and in Section [4] we discuss and comment 
our results. Finally we summarise our work and draw con¬ 
clusions in Section [21 In the Appendix El we present a per¬ 
turbative approach to describe the effect of non-local bias 
on third-order statistics in configuration space. 


2 SIMULATION 


Our analysis is based on the Grand Ghallenge run of the 
Marenostrum Institut de Giencies de I’Espai (MIGE) simula¬ 
tion suite to which we refer to as MIGE-GG in the following. 
Starting from small initial density fluctuations at redshift 
2 = 100, the formation of large scale cosmic structure was 
computed with 4096® gravitationally interacting collision¬ 
less particles in a 3072 h“^Mpc box using the GADGET - 2 
code llSpringelll2005ll with a softening length of 50 /i“^kpc. 
The initial velocities and particle displacements were gener¬ 
ated using the Zel’dovich approximation and a GAMB power 
spectrum with the power law index of ris = 0.95, wich as nor¬ 
malised to fulfil (Ts = 0.8 at z = 0.0. The cosmic expansion 
is described by the AGDM model for a flat universe with 
a mass density of Dm = ^dm -b Dj, = 0.25. The density of 
the baryonic mass is set to Dt = 0.044 and Q.dm is the dark 
matter density. The dimensionless Hubble parameter is set 
to h = 0.7. More det ails and valida t ion tes t on this simu¬ 
lation c an b e found in FosaIbaet_alJ l|2015al ), iFosalba et al.l 
ll2015bl ) and IHoffmann et al.l (l2015bli . 

Da rk matter haloes w ere identified as Friends-of-Friends 
groups dDavis et al.lllOSSl 'l with a redshift independent link¬ 
ing length of 0.2 in units of the mean particle separation. 
These halo catal ogs and the correspo nding validation checks 
are presented in iGrocce et al.l ll2013lf. _ 

As in our previous analysis ( Hoffmann et al.ll2015lJl we 
divide the haloes into the four redshift independent mass 
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Table 1. Halo mass samples. Np is the number of dark matter 
particles per halo, Nhalo is the number of haloes per sample in 
the comoving output at redshift z = 0.5. 


sample 

mass range [lO^^h ^Mq] 

Np 

^halo 

MO 

0.58 - 2.32 

20 - 80 

122300728 

Ml 

2.32 - 9.26 

80 - 316 

31765907 

M2 

9.26 - 100 

316 - 3416 

8505326 

M3 

> 100 

^ 3416 

280837 



12 13 14 15 

log( Mass [h“'l\/ls] ) 


Figure 1. Peak-background split prediction of the bias para¬ 
meters bi, C 2 and C 3 for the mass samples M0-M3 versus the 
mean halo mass of each sample at redshift z = 0.0 and z = 0.5 
(diamonds and triangles respectively). The mass ranges of the 
samples are marked by vertical dashed lines. 

samples MO, Ml, M2 and M3, specified in Table [T] These 
samples span a mass range from Milky Way like haloes up 
to massive galaxy clusters. For studying the bias estimators 
we use dark matter particles and haloes identified in the 
comoving outputs at redshift z = 0.0 and 0.5. 

2.1 Bias predictions 

The peak-background split model predicts that the 
bias parame ters in equati o n (fTjl are a function of the 
halo mass (iBar deen et ^ llOsd : I Cole and Kaiser! 1 19891 : 
ISheth and Tormenlll999l ). For deriving these predict i ons w e 
fit the mass function with the model of iTinker et al.l ll2010l) . 
The low mass sample Ml is thereby excluded from the fitting 
range, since we expect the mass function measurements to 
be a ffected by noise in the halo detection in this mass range 
(see iHoffmann et al.ll201^ . for a detailed discussion). How¬ 
ever, this exclusion does not affect the sign of the predicted 
bias coefficients. The mean bias coefficients of each mass 
sample, shown in Fig. is then obtained by weighting the 
bias prediction with the normalised halo mass distribution. 

In Fig. m we show the predicted m ass dependency 
for bias coefficients up to third order. In IHoffmann et al.l 


(l2015ah we compare the predictions for the linear bias para¬ 
meter with measurements from two-point cross-correlations 
and find an agreement at the 10% level. Comparing the 
second- and third-order bias predictions to measurements 
from third-order statistics and a direct analysis of matter 
and halo density fluctuations we find a good qualitative 
agreement (Bel, Hoffmann, Gaztanaga, in preparation). We 
shall use these predictions for discussing our results in Sec¬ 
tion [21 


3 BIAS ESTIMATORS 

In this section we study various bias estimators from second- 
and third-order clustering statistics of haloes and matter, in 
order to quantify and understand differences between these 
estimations. Such an understanding is crucial for using third- 
order statistics in order to break the degeneracy between the 
linear galaxy (or halo) bias and the linear growth of ma tter 
fluctuations, as we discussed in (IHoffmann et al.ll2015bl) . 

In this previous study we found that the linear bias from 
the reduced three-point correlation in configuration space 
tends to overestimate the linear bias from the two-point 
correlation, even when the analysis is performed at very 
large scales (> 30 /t~^Mpc). Similar findings have been re - 
ported in the literature (e.g. iManera and Gaztanagall201ll ). 
Such deviations can be ex pected from non- l ocal contribu¬ 
tions to the bias function dchan et al ] |2012l:lBald^ et al.l 
I2OI2I) . Furthermore, non-linear terms in the perturbative 
expansion of correlations functions, which are usually neg¬ 
lected in the analysis of clustering measurements, can con- 
tribut e to the deviations between the different bias estim¬ 
ators (| Pollack et al.l 1^121 ). The goal of this section is to in¬ 
vestigate the effect produced on various bias estimators of 
non-linear (sub-section [3T} and non-local (sub-sections 13.21 
and ITSI) contributions to the biasing function. 

Throughout this investigation we shall use the follow¬ 
ing notations: the density contrast of a stochastic field at 
position ri will be shortly written as Sy^i = (5;/(ri) where y 
refers to either haloes (h) or matter (m). Note that, when 
no confusion is possible, we will drop the index i referring 
to the position. 

3.1 Bias from second-order statistics 

We study second-order statistics of the halo- and matter 
density fields in terms of the two-point correlation and the 
variance. The density fields are therefore smoothed with a 
spherical Top-hat window function of radius R. The two- 
point correlations can be defined as the average products of 
density contrasts at the positions ri and r 2 , 

r" = (5.,i5,,2)(ri2), (4) 

which is a function of the distance ri 2 = |r 2 — ri |. If the two 
density contrasts belong to the same density fields (haloes 
or matter, i.e. x = y) equation (|3)) defines the two-point 
auto-correlation, denoted in the following as ^ and for 
the matter and the halo field respectively. Note that, for a 
given smoothing window, the variance corresponds to the 
two-point correlation in the limit where ri 2 = 0. In the case 
of two different density fields (halo and matter) equation (jd)! 
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defines the two-point cross-correlation, which we denote as 
in the following. 

The average (...) refers to all pairs, independently of 
their orientation. We thereby follow the common assump¬ 
tion of ergodicity such that the ensemble average can be 
estimated from the mean over the comoving volume in each 
simulation snapshot. The two-point correlation is therefore 
an isotropic quantity, in contrast to the three-point correla¬ 
tion, which we study later in Subsection 13.21 Note that we 
ignore redshift space distortions in this study. 

If the biasing relation between matter and haloes was 
linear {Sh ~ biSm) then the halo variance = (S^), the halo 
two-point auto-correlation = {Sh,iSh, 2 ), the halo-matter 
cross-variance = (ShSm) and the halo-matter two-point 
cross-correlation G = {Sh,iSm, 2 ), would provide equivalent 
estimators of the bias parameter b. Hence, we can define four 
linear bias estimators based on second-order statistics as 


and (gge 

iBernardeaul Il99a l. The non-linear terms in present in 
equation and (nsj, differ from each other. This means 
that non-linear bias affects the auto- and cross-correlation 
function of haloes differently. The corresponding expansions 
for the variance and the cross-variance of haloes can be 
straightforwardly obtained by taking the one-point limit in 
the equations (0) and (nni) (i.e. ri 2 = 0). Hence, the reduced 
correlator Cnk corresponds to the reduced cumulant Sn+k 
and the two-point correlation function (auto and cross) con¬ 
verges to the variance (auto and cross). These two expan¬ 
sions show that, at third order in cr^ and one needs to take 
into account bias parameters up to fifth order (cg). However, 
since we focus our analysis on scales at which the matter 
density contrast is small (Sm <C 1), we will approximate the 
bias function as a third-order polynomial. In the following, 
we therefore set the fourth- and fifth-order parameters to 
zero. 


bn 

= a-h/o-, 

( 5 ) 

bi 

_ 2/2 

( 6 ) 

h 

= 

( 7 ) 

lX 


( 8 ) 


Any deviation from the equivalence of these estimators 
provides a way to study the impact of possible non-linear 
terms in the biasing relation between matter and haloes (e.g. 
iManera et aLllioKl : ISaito et al.ll20lJ . respectively). In fact, 
if the biasing function is a local relation, defined by its Taylor 
expansion (see equation ([T} ), one can approximate the auto- 
and cross-correlation of haloes as a f unction of one- and 
two-p oint statistics of the matter f ield. iFrv and Gaztanagal 
(Il993tl and iBel and Marinonil (l2012f l provide the correspond¬ 
ing expressions respectively for one- and two-point cumu- 
lants up to fifth order (i.e. a-nd n + m ^ 5), keeping 

terms up to second order (leading order plus one) in terms 
of ^ and a^. 

Since we focus on non-linear effects, we express the 
second-order statistics of haloes (auto and cross) up to two 
orders after the leading order, which is keeping only terms 
of third-order in ^ and cr^ and lower. In the case of the two- 
point auto-correlation we obtain 


~ 5 -I- (C2C12 -f 03)^0-^ + 

+ U ^3+C4(— + -) + - 

+ -^C22 + y C3C12 -I- 

+ (C2C3C12 -I- y <= 4 ^ -I- ( 9 ) 

while we find for the cross-correlation 

^ + -( 02 ( 712 + 03 )^ 0 -^ 

where Cnk and S„ are respectively the reduced correlat¬ 
ors and cumulants of the matter field, defined as Cnk = 


The equations and GQl) show that, if only the lead¬ 
ing order terms are kept, then the bias estimators 
converge to the linear bias bi. However, they also tell that, 
if higher order terms are contributing to the signal, then 
significant deviations between the estimators must be de¬ 
tectable. In the following, we exploit this fact to test the 
impact of non-linearities on the estimators ®-®. 


The auto- and cross-correlation funct ions are estimated 
by ap plying the estimator proposed by iBel and Marinonil 
(l2012ll on the halo and the dark matter particle distribu¬ 
tions. Its implementation in the M ICE-GC simulation is de¬ 
scribed in lHoffmann et ^ (l2015lJ) . to estimate the errors we 
use a Jack-knife method on 64 cubical cells. For studying the 
scale dependence we vary the smoothing scale R, while we 
keep the correlation length ratio ri 2 /i? = 2 fixed. 


The bias measurements for the four mass samples MO¬ 
MS at the redshifts z = 0.0 and 2 = 0.5, derived from 
the estimators (IS])-® are presented in Fig. [2] We can see 
that, on the considered scale range, the estimators involving 
two-point statistics are in good agreement with each other 
and display negligible scale dependence. This indicates that 
higher-order terms in the equations m and GSJ can indeed 
be safely neglected. However, estimators involving the vari¬ 
ance differ significantly from b^ and bf , especially at small 
scales when the matter variance approaches unity. Fur¬ 
thermore, they do not mutually agree, even on large scales . 
Sim ilar results have been obtained bv iManera et ah] (l2010ll 
and iManera and Gaztanaga l^20 11^. In order to explain dif¬ 
ferences between b^ and b^ Manera et al. Ilm^) set C3 = 0 
and drop the term in into the higher order expansion; al¬ 
though they did not have a theoretical motivation for that. 


Following their approach and using the equations (0 
and (HSJ we expand the expressions © - ® up to the 
second-order in and ^ and express ba, b^ and with 
respect to bf 
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z = 0.0 z = 0.5 



Figure 2. Comparison between four linear bias estimations from 
measurements of the auto- and cross- variance and the two-point 
auto- and cross- correlation (blue, green, red and black respect¬ 
ively), derived via equation Q - ([8jl. Results for the redshifts 
2 : = 0.0 and z = 0.5 are shown in the left and right panel respect¬ 
ively. The different smoothing scales R of the halo and matter 
density fields, used for the bias estimations, are shown on the 
a:-axis. For measuring ^ and we fix r \2 = 2i?. 



( 


bi 


2 



( 


bi 


2 



C2(<S'3 — C12) + 


+ 


j(S4 - C?2) + j(S4- C13) 


+ ^(353-Ci2) + | 


C2(<5'3 — Cl2)cr^ + 


+ ^(S4-Cl3) 


+ 


C2C3 


(5'3 - C12)] 


|?+|[C22-C?2] 


( 11 ) 




( 12 ) 


+C2C3C'l2Ccr^ + 

6 


(13) 


With the expressions above one can partially explain the 
deviations between the different bias estimators, displayed 
in Fig. O We thereby rely on the fact that differences 
between the linear bias and differences between the second 
power of these values are monotonically related to each- 
other, due to the positive amplitude of the linear bias. The 
first interesting point is to notice that equation m does 
not exhibit any contribution from at the leading or¬ 
der. Moreover the contribution from vanishes since it 


has been show n that correlators obey the factorisation rule 
Cnk = CinCik teernardeaull 199di . As a result equation (fT3t 
contains only terms in which is small on the considered 
scale range (ri 2 > 20/i“^Mpc). This explains the good agree¬ 
ment between b^ and b^ for all halo mass ranges. 

The exp ression which we obt ained for 65 differs from the 
one given bv iManera et al.l (l201Clf) presumably due to a trun¬ 
cation of the higher-order expansion of the halo two-point 
cross-correlation function by these authors. In addition, our 
results seem to confirm the necessity they had to set the 
contribution of order cr'^ and the third-order bias C 3 to zero 
in order to explain their results. Note also that at very small 
scales, Fig.[2]shows that b^ tends to be larger than b^. This 
small tendency can be explained by the fact that the leading 
order term of equation m is positive. 

Regarding the comparison between b^ and b ^, the lead¬ 
ing order term in expression m shows that their relative 
amplitude depends on the sign of C 2 . If it is negative b^ is 
lower than b ^, which is the case for the mass samples MO- 
M2 at 2 ; = 0 and MO-Ml at 2 = 0.5. The peak-background 
split bias predictions, displayed in Fig. [T] shows that C 2 is 
indeed expected to be negative for these mass samples (see 
Section [ 2 ] for details about the prediction). The deviations 
between the bias from the variance and the two-point cross¬ 
correlation is, at low masses, in good agreement with the one 
predicted by equation m- In fact, the leading order term 
contains two contributions, one in and another in C 2 . The 
latter being dominant for low masses where |c 2 | < 1. It fol¬ 
lows that ba is expected to lie between and b^ which 
is confirmed by Fig. [2] However, for higher mass bins we 
observe an opposite tendency to the expected one. By com¬ 
bining equations uni and m one obtains that the relative 
position between b^ and b^ is 



C2 2 , 


C2C353 + f 

4 6 


(14) 


As a result when C 2 > 0, the bias from the variance is ex¬ 
pected to be greater than the bias from the cross-variance. 
Note also that, even in the particular case of C 2 — 0 (M2 
at 2 = 0.5), the remaining term in cl shows that foo- must 
be greater than The measurements, shown in Fig.[2l ex¬ 
hibit a different behaviour than the one just described. In 
fact, for the high mass bins M3 at 2 = 0 and M2, M3 at 
2 = 0.5, the estimate from ba stays below the one from foj. 
The bias from the variance apparently saturates to the value 
of the bias from two-point correlations and therefore seems 
to constitute a more robust linear bias estimator than the 
cross-variance. Note that, this tendency is observed even at 
very large scale R = 40 /i“^Mpc which suggests that non¬ 
linear corrections are not responsible. 

We conclude that higher-order contributions arising 
from non-linear bias can only explain partially the relative 
amplitudes between the linear bias estimators displayed in 
Fig. [21 This result motivates the introduction of a non-local 
componen t in the biasing relation, as suggested in the liter- 
ature le.g. lChan et al.ll 2012 l : [^aldauf et al.il 201 j : ISaito et al.l 
l2014h and also shows that one needs to go beyond second- 
order statistics. However, we do not expect non-local bias 
to significantly affect our second-order halo statistics, espe¬ 
cially because of the large smoothing scales applied in this 
analysis. This might not be the case for third-order statist- 
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ics. We therefore speculated in iHoffmann et ^ (l2015lJ) that 
a non-local component in the biasing relation could explain 
the differences between second- and third-order bias estim¬ 
ators. In the next subsection, we therefore generalise our 
analysis in studying linear and quadratics bias estimation 
from third-order statistics (auto and cross). Finally, we note 
that the bias estimator based the two-point cross-correlation 
function {bf ) can be considered as a reliable estimator of the 
linear bias for, we thus consider it as a reference in the fol¬ 
lowing. 


3.2 Bias from three-point cross-correlation bg 

The three-point correlation of the three generic density fields 
that we label x, y and 2 at the positions ri, r 2 and rs 
(respectively) is defined as 


= (5.,i5,.2<5.,3), (15) 

where the vectors ri, r 2 and rs form a closed triangle which 
can be parametrised in terms of the size of its three legs 
Xij = |rj — ri| or in terms of the two legs ri 2 , ris and the 
angle 023 = acos(fi 2 • fis) between them. As for the two- 
point correlation in Section [33](---) denotes the mean over 
all triangles in the analysed volume. However, in contrast to 
the two-point correlation, where the average is taken over all 
pairs of 5 independently of their direction, C, is not isotropic 
as it is sensitive to the shape of the large-scale structure. It 
therefore provides access to additional information. In case 
of the auto-correlation the three density contrasts in equa¬ 
tion m refer to the same density field. The correspond¬ 
ing three-point correlation function of haloes and matter 
are therefore and respectively. To measure the 

bias from the three-point cross-correlation between haloes 
and matter densities we compute ^ ^mhm^ ^mmh^ These 
quantities are the n compared to the hierarchical three-point 
cross-correlation (lFrvlll98^ f 


j-hm _ ^hm ^hm . ^mh ^hm . ^mh ^mh 

SH =^12 S13 + S12 ^23 + S13 S23 • 


(16) 


Note that here referrers to the two-point cross¬ 

correlation between haloes at position ri and matter at 
position rj, which is called G in fh® remainder of this 
article. Combining equation GSl) and GSl) one can define 
the s ymmetric reduced three-point cross-correlation func¬ 
tion ( Pollack et al.|[^12ll 


1 /-hmm I /-mhm i /-mmh 

— cr • 

The reduced three-point auto-correlations for matter and 
halo density fields are defined analogously as Qm = 
and Qh = • This way Qh and 

quanti fy any departure from the hierarchical ansatz (IFrvI 
Il984bl ). In the following we will refer to the reduced three- 
point correlation as the three-point correlation. 

Inserting the non-local quadratic bias model (equation 
([^) into the definition of the three-point correlation for ha¬ 
loes yields, via a second-order perturbative expansion, in the 
limit of small density fluctuations and large triangles 



« Qm 


Figure 3. Left, three-point auto- and cross-correlation for mat¬ 
ter and haloes (Q/i, and Qm) in the mass sample M2 at 
redshift z = 0.5, measured using triangles with fixed legs of 32 
and 64 /i“^Mpc for different opening angles a (circles, crosses 
and squares respectively). Fits based on the local bias model (i.e. 
g 2 = 0) to Qh and from equation CEJ and esj are shown as 
thick dashed-dotted and dash-double dotted lines. Predictions for 
the matter three-point correlation and the non-local component 
from perturbation theory {Q^ and shown as thick 

and thin dashed lines respectively. Right the and versus 
Qm relation, used for deriving the linear and quadratic paramet¬ 
ers bi and C 2 in the local bias model. Thin dashed-dotted and 
dash-double dotted lines show the non-local contributions to Qh 
and Q^ respectively, using g 2 from AQcg, equation 


Qh — ~g~ {Qm + [C2 + y2Q'nloc\ } , (1^) 

which can be generalised to the case of three-point cross¬ 
correlation, 


Qh — T + g2Qnioc]} (19) 

(see Appendixl^ . These expressions differ signihcantly from 
the ones obtained from the local bias model, as they include 
the non-local contribution to the three-point halo correlation 
Qnioc, which we present in more detail belo w. The expression 
for the local model, which we assumed in iHoffmann et al.l 
l|2015bh corresponds to a vanishing non-local bias para¬ 
meter, i.e. g 2 = 0. Since Qnioc is a function of the opening 
angle a = 0 : 23 , the three-point halo auto- and halo-matter 
cross-correlations are therefore no longer linearly related to 
the matter three-point correlation. This a dependence arises 
from the fact that Qnioc originates from tidal forces, which 
modify the shape of matter fluctuations. 

We can predict Qnioc from the power spectrum, assum¬ 
ing that the perturbations of the density held 5m are small 
(i.e. the divergence of the velocity held Ov in equation m is 
linear and therefore equal to 5m- It also requires to assume 
that the legs of the considered triangles are large compared 
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Figure 4. Linear bias of the halo mass sample M2 at the redshift 
2 = 0.5, measured independently of second-order contributions 
(local and non-local) to the three-point correlation via equation 
(123 by combining and Qm (op^n circles). The fit to 

the measurements is shown as solid line. The linear bias meas¬ 
urements, derived from Qh and within the local bias 
model (via equation (Hi, lO and (Hii respectively with g 2 = 0) 
are shown as dashed, dash-dotted and dash-double dotted lines 
respectively. 


to the smoothing radius R (large separation limit, i.e. in 
equation Q, W{Rqi 2 ) W{Rqi)W{Rq 2 ))■ For such con¬ 
ditions the perturbation theory (hereafter also referred to 
as PT) offers the possibility to express the non-local three- 
point correlation Qnioc in terms of 


r 123 


^(^ 12 ) + 3 


<(>'(^12) 

ri 2 


where 


C(ri3) -I- 3 


ri3 


L2 (cos 012), 


( 20 ) 


(t){r) = 






kr 


and <(>^(r) = ^{r). One can show that 


( 21 ) 


^ 2 f ri23 + r3i2-I-r23i 

Qnloc{'r'l2,ri3,a) ^ g 

The angular dependence of the non-local component of the 
three-point halo auto- and cross-correlation functions is en¬ 
coded in equation (I20II via the seco nd-order Legendre poly¬ 
nomia l L 2 (cosqi 2 ). As shown by iBarriga and Gaztanaaa 
l|2002h in their equation (8), at the tree-level and for large 
separations, the matter three-point correlation can be ex¬ 
pressed in the same way as expression That is with re¬ 
spect to circular permutations of a function ri 23 , expressed 
as a monopole, a dipole and a quadrupole in cos(qi 2 ) (sim- 
ilar Legendre expansio n have been used in Fourier space by 
ISchmittfull et aDl201!Th As a result, a non-local component, 
such as the tidal field Q 2 (equation ([S])), modifies the amp¬ 
litude of the quadrupole of F 123 by an amount proportional 
to the non-local bias g 2 (see Appendix 

Moreover, by comparing equations nil) and (HU) one 


can see that quadratic and non-local contributions to the 
halo three-point correlation affect the cross-correlations by 
a factor 1/3 less than the auto-correlation. The linear bias, in 
contrast, affects the auto- and cross-correlation equally. We 
will use this property to isolate the linear from the quadratic 
and non-local bias, as explained below. 

In the following we study non-local contributions to 
the halo bias in the MICE-GC simulation. Ours measure¬ 
ments of the three-po int correlation are base d on t he al- 
gori thm suggested in ISarriga and Gaztahagal (l 2002 l ~l (see 
also iHoffmann et al.ll2015bl . for studies of numerical effects 
in this algorithm and the impact of covariance between an¬ 
gular bins on the bias measurements). Errors are derived 
from 64 cubical Jack-Knife samples. We first focus on the 
mass bin M2 at redshift z = 0.5 to present our methods 
for extracting the parameters 61 , C 2 and g 2 from three-point 
correlations, which were computed using triangles with fixed 
legs of ri 2 = 32 and ri 3 = 64 h“^Mpc. Afterwards we will 
present results for all mass samples and redshifts at various 
scales. 


3.2.1 Local bias 

Our first method for measuring bias from three-point cor¬ 
relations is based on the local bias model (52 = 0). The 
linear and quadratic bias parameters are computed from 
the equations (HU and m by fitting the 61 and C2 para¬ 
meters which allows for mapping Qm into Qh and Q^, i.e. 
Qh = (Qm + C 2 )/bi an dQ/ = (Qm + ^)/5i. T he fitting pro¬ 
cedure is explained in IHoffmann et ^ (|2015lJ l. The two es¬ 
timations of the doublet 61 and C 2 are respectively called 6 q, 
cq and , Cq. In Fig. |3]we show how well a linear relation, 
expected from the local bias model, describes the mapping 
between the matter three-point correlation and the three- 
point auto- and cross-correlation functions of haloes. How¬ 
ever, we can see that the slope of the linear relation is differ¬ 
ent when considering auto- and cross-correlations, which in¬ 
dicates that the two methods deliver inco nsistent resu l ts (see 
right panel of Fig. 0 . As explained in iGhan et al.l (I 2 OI 2 II 
this linear relation between matter and haloes might arise 
from a projection effect due to the fact that we neglect the 
non-local component g 2 Qnioc- In Fig. [3]we show that, if the 
contribution of Qnioc is small compared to Qm (i.e. 32 is 
small, see Fig. B, then they can indeed be approximately 
related by a linear relation. The ignored non-local contri¬ 
bution to halo three-point correlations can therefore be ab¬ 
sorbed by 61 and C2, without substantially decreasing the 
goodness of the Qh and Q^ fits. The same eff ect has been 
shown in Fourier space bv iBaldauf et ah 1 (1^) in their Fig. 
1. This ignorance might lead to incorrect bias measurements, 
unless g 2 — 0. 


3.2.2 Non-local bias 

Our second method for measuring bias from three-point 
correlations is a new a approach, which combines auto- 
and cross-correlations. These two statistics can be combined 
in two different ways which allow us to isolate the linear 
bias from quadratic and non-local contributions to the bias 
model. Both combinations take advantage of the fact that 
the linear bias, 61 , affects Qh and Q^ equally, which is not 
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^ ^ nloc 

Figure 5. Left: The difference between the three-point auto- and 
cross-correlation and (AQcg), multiplied with the linear 
bias from the two-point correlation (open circles) as a function of 
triangle opening angle a for M2 at z = 0.5. According to equa¬ 
tion ilM this expression is equivalent to C 2 + 92 Qnloc- If bias 
function is quadratic and local, the measurements should corres¬ 
pond to the quadratic bias parameter C 2 , as measured from the 
halo-matter cross-correlation via equation m with g 2 = 0 (dash- 
dotted line). Fits from equation |[24t . based on the PT prediction 
for the non-local component, Qnioci are shown as thick solid line. 
The quadratic bias parameter, derived from this fit (AQcg) is 
shown as dashed line. The corresponding effective quadratic bias 
parameter (equation OB) is shown as dashed-double dotted line. 
The three-point correlations were calculated from triangles with 
fixed leg of ri 2 = 'ri 3/2 = 32 h“^Mpc. Right: same as left panel, 
but showing the measurements versus the Qnioc prediction for 
each opening angle. 

the case for the quadratic and non-local contributions, C 2 
and 92 . The linear bias can be obtained by combining the 
equations m and dm, 

bAQ = -2||, (23) 

where AQ = Qh — 3Q^. The interesting property of this 
linear bias estimator is that it is independent from quad¬ 
ratic (local and non-local) contributions to the bias function. 
It can therefore be used to verify if such contributions are 
indeed the reason for deviations between linear bias estim¬ 
ations fro m two- and three-point correlations, as we spec¬ 
ulated in iHoffmann et ^ (l2015bli . Note that the relevant 
quantities involved in equation (1231) depend on the open¬ 
ing angle a, so does the estimator h^Q- Hence, our final 
&AQ measurement is derived by htting a constant to hAqia)- 
The use of AQ has also the advantage that off-diagonal ele¬ 
ments in the covariance matrix between different opening 
angles are smaller. This covariance is difficult to access and 
its Jack-Knif e estimation can affect the bias estimation at 
the 5% level (IHoffmann et al.ll2015bl '). 

We show &aq for M2 at « = 0.5 in Fig. U at each angle 
probed by the three-point functions together with the cor¬ 
responding fit. In the same figure we also display the es¬ 



^ nloc 

Figure 6. Same as right panel of Fig. [5] but for the halo mass 
sample M3 at the redshift 2 = 0.0. Results are shown for triangles 
with fixed legs ri2 = ^13/2 of 12 and 36 /i~^Mpc (open circles 
and filled triangles respectively). Dashed and dashed-dotted lines 
are the corresponding fits to equation dsi, used to derive the 
quadratic and non-local bias C 2 and g 2 . The black solid line cor¬ 
responds to the local Lagrangian prediction for C 2 + 52Qn/oc with 
C2 = 0 and g 2 = 272/51, while 72 = — ( 2 / 7 )( 6 i — 1 ). The 72 values 
from these fits are shown in Fig. [ 8 ] 


timations for the linear bias, derived from three-point auto- 
and cross-correlations [hq and 5 qx , obtained from equations 
ca and GSl), assuming the local bias model, i.e. 32 = 0 ). 
As a reference, we also include the linear bias measurements 
from the two-point cross-correlation, , which we consider 
to be a reliable estimate of the true linear bias (see Section 
13.111 . The comparison in Fig. |4] re veals that the measurement 
and the fit of the hybrid bias estimator b/^q are consistent 
with the reference fej, while we see that the biases obtained 
from bq and 6 qx are over estima ting the linear b i as. Thi s 
result confirms our speculation in IHoffmann et ^ (|2015lJ l. 
that differences between 65 and bq are mainly due to a non¬ 
local term in the bias model. Furthermore we verified that 
also the magnitude of the overestimation is in agreement 
with what we expect from neglecting non-local contribu¬ 
tions. However, as we pointed out in this aforementioned 
study, linear bias measurements from three-point correla¬ 
tions can also be affected by numerical effects, originating 
from the algorithm employed for deriving Q as well as short¬ 
comings in the estimation of the Q covariance between dif¬ 
ferent angular bins. 

One can notice that, as expected in case of non-local 
bias, the overestimation is larger in case of the auto- 
esti mator {bq) compared to the cross-estimator ( 6 qx ). Note 
that IPollack et al.l (I 2 OI 2 I I found an opposite trend, analys¬ 
ing a different simulation in Fourier space with a different 
mass resolution and cosmology. Their linear bias measure¬ 
ments from the bispectrum are closer to peak-background 
split predictions than measurements from three-point cross¬ 
correlations, while the p redictions might be low er than the 
true linear bias (see e.g. IHoffmann et al.ll2015al l. 
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Figure 7. Scale dependence of bias measurements from three-point correlations. Measurements were derived using triangles with 
^' 13/^12 = 2 configuration, while scale of the smaller triangle leg, ri 2 , is denoted on the x-axis (slightly shifted for clarity). Results 
are shown for the four mass samples M0-M3 (from left to right) at redshift 2 : = 0.5. The linear and quadratic bias (top and central 
panel respectively), obtained from three-point auto- and cross-correlations using the local bias model {g 2 = 0) via the equations (USD 
and 03 are shown as blue circles and green crosses respectively. The linear, quadratic and non-local bias, measured using combinations 
of three-point auto- and cross-correlations (AQ, equation (1231 ) and (I24D ) are shown as red triangles. The linear bias is compared with 
reference measurements from the two-point cross-correlation (5^), shown as black, dashed line. The quadratic bias measurements are 
compared with effective quadratic bias (equation 0371 1 h shown as orange triangles. The non-local bias parameter p 2 = 272/61 (bottom 
panel) is compared with predictions from the local Lagrangian model (72 = —(2/7)(6i — 1)). (red dashed line). Error bars denote Icr 
uncertainties. 

3.2.3 Linear and quadratic terms 

In order to further verify the presence of non-local contri¬ 
butions to the bias function we separate Qnioc from Qm by 
subtracting from Qh and define 

AQcg = Qh ~ Qh ~ Q^qQ^iIoc] 5 (2^^) 

where caq and qaq are the estimators of C 2 and g 2 in equa¬ 
tions dm find . If tliG bifis function is c^uudrfLtic in Sj^i 
and local then the non-local bias parameter g 2 is zero. Hence, 


CAQ = biAQ,g{3/2) (25) 

should correspond to C2, independently of the opening angle 
a . In the left panel of Fig. [ 5 ] we show caq, together with 
C2, estimated from Qh and (cq and Cqx respectively) 
from the local bias model (equations (1181) and (1191) . for g 2 = 
0). In the same figure we also show C 2 derived from fitting 
AQcg taking into account a possible non-local bias g 2 - The 
first important point is that the measured caq shows a very 
clear angular dependence, with a maximum at around 80 
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degrees and a positive curvature. The local quadratic model 
therefore clearly fails in describing the impact of bias on 
three-point correlations. 

The right panel of Fig. [5] displays the measurements of 
AQcg with respect to the tidal, non-local component of the 
three-point function Qnioc, as predicted by equation (El). 
It shows that measurements are compatible with a linear 
relation between AQcg and Qnioc, which is expected in the 
second order non-local bias model. It therefore confirms the 
presence of non-local bias due to tidal forces {G 2 ) in the mat¬ 
ter field. For such a linear relation, the value of AQcg corres¬ 
ponds to the second-order bias caq when Qnioc = 0 , while 
its slope provides a direct insight to the non-local bias g 2 - 
Finally, by comparing the various estimates of the second- 
order local bias C2, one can see that, when assuming a local 
bias, the results differ significantly from the ones obtained 
when taking into account a possible non-local component 
into the biasing relation. However, the estimation of caq 
and g 2 might be affected by shortcomings of the PT predic¬ 
tions for Qnioc- 

To verify how well the PT predictions describe the 
measurements at different scales we show the relation 
between AQcg and Qnioc derived from triangles with fixed 
legs ri 2 = r\zl2 of 12 and 36 h~^ Mpc respectively in Fig. 
for the mass samples M3 at z = 0.0, for which we find a large 
non-local bias amplitude. At large scales (ri 2 = 36 h~^ Mpc) 
the slope of the measured AQcg - Qnioc relation is compar¬ 
able with the local Lagrangian prediction. Interestingly at 
small scales the AQcg - Qnioc relation is also linear, while 
the slope has the opposite sign than at large scales. The 
linearity at small scales indicates that higher-order terms 
enter the Q in a similar way as second-order non-local con¬ 
tributions to the bias function. This suggests that linear 
bias measurements can be improved by using the prediction 
for Qnioc, while using the local Lagrangian prediction for 
the non-local bias g 2 is only appropriate at extremely large 
scales. 


3 . 2.4 Scale dependence 

Based on the methods for measuring linear, second-order 
and non-local bias from Q and Q^ {Bq, Bqx , Baq, cq, Cqx 
and Caq), which were presented above, we now apply our 
analysis to each of the mass samples M0-M3 at z = 0.5. 
We study the scale dependence of our results as before by 
varying the size of the triangle leg ri 2 between 12 /i“^Mpc 
and 36 h“^Mpc) while fixing ri 2 /ri 3 = 1/2. The various 
bias estimations are presented in Fig.[7]for different triangles 
sizes, defined by the length of ri 2 . 

From the comparison between the different bias estima¬ 
tions we draw similar conclusions as from the example of M2 
at z = 0.5. On linear scales (sufficiently large triangles) the 
linear bias parameters obtained from each method reaches 
a regime in which they become scale independent. However, 
they do not converge to the same value. In case of the linear 
bias, only Baq is in agreement with hf , while bq and bqx 
overestimate the linear bias; this overestimation is stronger 
in the case of Q than Q ^. 

The scale dependence shown for the high mass bin M3 
in Fig. [7] shows that, if the analysis is performed at too 
small scales, then one can measure a positive non-local bias, 
while it is in reality negative (as we see from the results 


derived at large scales). We indeed verified that for highly 
biased tracers and small triangles, the curvature of caq flips 
from positive to negative. This scale dependence indicates a 
domination of the signal by non-linear terms in these cases. 
However, for lower mass samples and large scales, for which 
we expect non-linear contributions to converge to zero, we 
still see a strong angular dependence of caq, which speaks 
for the presence of non-local bias contributions. Hence, we 
have shown that all the mass sample used in the present 
analysis exhibit a detectable non-local component. 

In case of estimators of the second-order bias C 2 (middle 
panel of Fig. 0 we observe the same tendency, however we 
do not have a reference estimate. As a result we shall be 
more confident in the estimation of C2 coming from caq. 
We compare the latter to measurements from third-order 
moments in Section [T] A comparison with C 2 derived from 
various methods will be presented in Bel et al. (in prepara¬ 
tion). 

The bottom panel shows that each mass sample com¬ 
prises non-local bias which significantly differs from the local 
model 72 = 0. These measurements therefore constitute the 
first detection of non-local bias in configuration space. In 
the case of MO and Ml the amplitude of the non-local bias 
strongly differs from the local Lag rangian biasing relat ion 
between the halo and matter field (IMo and Whitef[l996f l. 


3.2.5 Non-local to linear bias relation 


Following IChan et al.l ll2012ll . we compare our non-local bias 
measurements to the linear bias derived from the two-point 
cross-correlation in the top panels of Fig. [ 8 ] This compar¬ 
ison includes measurements at redshift 0.0 and 0.5 which 
are based on triangles with r\ 2 /r\s = 1/2 configurations. 
For very large triangles (ri 2 = 36 /i“^Mpc) our results in¬ 
dicate a linear relation between the non-local 72 and the 
linear 61 bias, as expected for the local Lagrangian bias¬ 
ing. However, the amplitude of this relation lies below the 
local Lagrangia n prediction, which is the opposite of what 
was reported bv IChan et al.l ll2012ll . Some work is currently 
ongoing, aiming to explain whether t hese differences result 
from the fact that IChan et al.l (l2012li conduct their meas¬ 
urements using the Bispectrum in Fourier space, while we 
employ the reduced three-point correlation in configuration 
space. A further contribution to the discrepancies could arise 
from differences in the simulation, such as mass resolution 
effects, or differences between cosmological parameters. 

As a matter of fact, our measur ed hi — 72 r elatio n 
shows the same ten dency as those o f Sheth et al.l ll2013ll . 
iBaldauf et al.l (l 2012 l l and ISaito et al.l ll2014h who also find 
(in Fourier space) 72 to be below the local Lagrangian pre¬ 
diction. The simulations employed in the two latter stud- 
ies are based on cosmologies similar to MICE-GC, while 
IChan et al.l (l2012h study a simulation with a initial power 
spectrum which is significantly different in terms of its nor¬ 
malisation 0 - 8(2 = 0.0) and spectral index Us. Note also 
that the departures from the local Lagrangian prediction in 
Fig.0are strongly scale dependent for highly biased samples 
(foj > 2 ), which indicates the p resence of non-lin ear contam¬ 
ination to Qm and Qnioc (e.g. ISaito et al.|[2014l L 
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Figure 8. Non-local bias parameters 72 and g 2 = 272/^1 versus 
the linear bias from the two-point cross-correlation, . Meas¬ 
urements, derived from AQcg = {Qh ~ Q^) via equation I I24II 
are shown as symbols with Icr errors. Results are shown for Q 
measurements based on triangles with fixed legs of ri 2 = 12, 24 
and 36 /i~^Mpc and ri 2 = ^ 13/2 configurations (from the top 
to the bottom). Crosses, open squares, open circles and open tri¬ 
angles show the non-local bias, measured for the mass samples 
MO, Ml, M2 and M3 respectively at redshift 2 : = 0.0. Stars, 
closed squares, closed circles and closed triangles show the corres¬ 
ponding M0-M3 measurements at 2 = 0.5. The measurements for 
ri 2 = 36 /i~^Mpc are approximated with 72 = —(2/7)(5i — 0.8) 
(black dash-dotted lines) and compared with the local Lagrangian 
prediction (- ^2 = —(2/7)(6i — 1)) as well as to a fit given by 
IChan et al. I mi± 72 = ~(2/7)(f)i — 1.43)), shown as solid and 
dashed black lines respectively. 


3.3 Bias from third-order moments 6^ 

An alternative approach for exploring the higher-order stat¬ 
istical properties of the large-scale structure are one- and 
two-point third-order statistics, respectively referred to as 
skewness and reduced correlator. They can be understood 
as the one- and two-point limits of the reduce d three-point 
correlation fe.g. lOaztanaga and FriemanlllQ!^ . Since these 
statistics are isotropic, they do not provide access to the full 
third-order hierarchy probed by the three-point correlation 
(i.e. the shape of the large-scale structure). However, this 
section will show that this apparent disadvantage leads to a 
cancellation of non-local bias which is useful for linear bias 
measurements, as we discuss later. 

The auto-skewness and the reduced auto-correlator of 
the matter density fluctuations 5m are respectively defined 
as 


and 



(26) 


Ci 2 (r) = 


{5m,l5m,2) 




(27) 


fsee lGoroff et ahlllQSd : [Bernarde'aulll996h . As in the previ¬ 
ous subsection a and ^ refer to the matter field. The auto¬ 
skewness and the reduced auto-correlator for halo density 
fluctuations, 5h, are defined analogously to equation (I26II 
and m- The skewness is directly related to the asymmetry 
of the one-point probability distribution of Sm, while the 
correlator tells how the quadratic field is correlated with 
itself on a given scale r. Their propert ies have been ex¬ 
tensively investigated in literature (e.g. BernardearJ 119961 : 
iGaztanaga et alT 20021 : ISel and Marinonif 20121) . As for the 
two-point and three-point correlation we also study the 
cross-skewness and the cross-correlator h order 
to investigate the impact of non-linearities, non-local bias 
as well as shot-noise on the measurements. We define these 
quantities as 


and 


c,x _ (ShSm) 

^ 3 ,h = 4 

^ X 


(28) 


'-'12,h = 


{5h,l5h,25m,2) 


(29) 


As in previous sub-sections the indices h and m refer, re¬ 
spectively, to the halo and matter density contrast. 

For measuring t he auto skewness S 3 and the auto cor¬ 
relator C 12 we follow ISel and Marinonil (l2012h by setting up 
a regular grid of spherical cells of radius R and counting the 
number of objects (haloes or dark matter particles) per cell. 
After assigning a number density contrast 5n to each grid 
cell we derive the auto-skewness S 3 as 


S3.JV 


(5%) - 3{5%)N-^ + 2N-^ 

{{5%)-N-^y 


(30) 


In order to estimate the reduced correlator C 12 , we consider 
the density contrast 5n,i = 5n for each grid cell. We then 
place an isotropic distribution of cells at distance r around 
the central cell and assign the number density contrast 5 n ,2 
to each of these surrounding cells. By averaging over all grid 
points in the simulation volume we estimate the correlator 
as 


{5n,i5% 2} ~ ‘2{Sn,i5n,2)N ^ 
(<5iv,i5iv,2>((5^) -V-i) 


(31) 


where N is the average number of halo or dark-matter 
particles per cell in the simulation box. These estimators 
are corrected for sho t-noise, assum ing the local Poisson pro¬ 
cess approximation llLavzerlll956h . Note that, in order to be 
able to handle the large number of dark matter particles, 
we use only 1/700 of the total number of particles in the 
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z = 0.0 z = 0.5 z = 0.0 z = 0.5 
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Figure 9. Left: reduced cross Skewness 5^, measured at the two redshifts 2 = 0.0 and z = 0.5 for matter (filled circles), and the four 
halo mass samples M0-M3 (diamonds, triangles, squares and stars respectively). Right: correspondi ng measurement s for th e reduced 
cross-correlator Ci 2 ‘ On each panel we display the tree-level PT predictions for 5 *3, equation (27) in | Hpffmann et ahl ll2015tJ) . For C 12 
we show two predictions. The dotted line displays the one from equation (28) in iHoffmann et "Si l|2015bl 'l and the solid black line shows 
the one taking into account the small separation contribution (Bel in preparation). We also display in grey dashed lines the fits to 
and C ^2 froiTi the equations l|35} and ilMli . 


dark matter simulation output. In principle this neglection 
of particles introduces additional shot-noise errors but we 
have tested that it does not affect the measurements. 

Following the same method, the cross-skewness is 
estimated as 


cx _ 

^3,h — 


{5N,hS%) — {SN,h5N)N 
i{SN,hSN)^ 


(32) 


where the sub-index h again refers to halo density field. 
Note that we correct the cross-skewness differently for shot- 
noise than the auto-skewness, since the former depends on 
the matter-field at the power of 2. Finally, since the cross¬ 
correlator does not involve any powers of the same density 
field larger than unity, the reduced cross-correlator h 
expected to be insensitive to Poisson shot-noise. Hence we 
estimate it directly from equation (1291) . The corresponding 
errors are estimated with Jack-Knife sampling, using 64 cu¬ 
bical cells. 

Our measurements of the skewness and the reduced cor¬ 
relator are presented as a function of smoothing scale R in 
Fig. [5]for the redshifts z = 0.0 and 2 : = 0.5. As in Sect ion IHTT] 
the scale dependence is explored by varying the smoothing 
scale R while we fix the ratio r/R — 2 in case of the cor¬ 
relator. In the same figure we also show the PT predictions 
for the dark matter field. The skewness of the dark matter 
field appears to be redshift independent in the linear regime 
[R > 20 / 1 “ ^Mpc), w hich is in good agr eement with the tree- 
level PT prediction (lBernardeaulll992 i. At smaller smooth¬ 
ing scales the redshift de pendence becomes more s i gnific ant 
which has been shown bv iFosalba and Gaztana^ 1 1998li to 
be due to non-linear loop corrections. The reduced correlator 
of the dark matter field shows only a weak redshift depend¬ 


ence, even for small sm oothing radii (R < lO/i ^Mpc). In 
I Hoffmann et al.l (l2015bli we have shown that the measured 
matter correlator is in disagreement with the tree-level PT 
prediction (dotted blue line in the right panel of Fig.|9}. We 
argued that this is due to the fac t that the standard tree - 
level prediction llBernardeaulll99^ : iBel and Marinonill201^ 
neglects the small separation contributions. We are indeed 
looking at the correlation between adjacent spheres. In the 
right panel of Fig.[9]we display also the tree-level prediction 
(magenta dotted line) which takes into account the small 
separation contribution (Bel in preparation) arising from the 
mode-coupling induced by the smoothing process. The ob¬ 
served agreement confirms that for small separations it is 
mandatory to take into account the smoothing in a proper 
way. 

In Fig. 1^ we also show for the first time meas¬ 
urements of the halo-matter cross-skewness and cross¬ 
correlator. Comparing these measurements to our results 
for the auto-skewness an d the auto-correlator, presented in 
I Hoffmann et al.l (l2015bli . we notice two qualitative differ¬ 
ences. The first difference is that the cross-skewness for the 
highest mass sample M3 increases at small scales, while we 
found an opposite trend only for M3 in case of the auto¬ 
skewness. The latter measurement should be more affected 
by shot-noise than the corresponding cross-skewness one or 
than measurements at lower mass samples; this is due to the 
low number density of M3 haloes. Furthermore the impact of 
shot-noise increases at small scales, where the cell volumes 
are small as well. The decrease of the M3 auto-skewness 
at small scales might therefore be attributed to a failure 
of Poisson shot-noise correction, indicating the presence of 
non-Poisson shot noise for highly biased samples. 

The second difference is that the amplitudes of the 
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cross-skewness and the cross-correlator increase monoton- 
ically with the halo mass. We did not find such a monotonic 
behaviour for the auto- skewness and the auto-correlator in 
Hoflhnami_etah I (l2015iJl which has been also confirmed by 


Angulo et al.l I 2008h . This difference can be explained by 


the fact that the bias parameters affect the auto- and cross¬ 
statistics differently, as we demonstrate in the following. 


3.3.1 Bias relations 

The skewness and the correlator of the matter held can be 
related to the corresponding halo and halo-matter cross¬ 
statistics via the non-local bias model. In analogy to the 
three-point correlation, one can express, at the leading or¬ 
der, the cross-statistics between haloes and matter as func¬ 
tion of the bias parameters by inserting the bias functions 
from equation © into the dehnitions (I26II - (I29II . In order 
to take into account the non-local component, we further 
assume that perturbations are well described by tree-level 
PT. It leads to the following expressions for the skewness, 
correlator, cross-skewness and cross-correlator of haloes 


typically above the measurements. This latter discrepancy 
is lower than in t he case of the auto-statistics, presented in 
I Hoffmann et al.l (l2015bll . We therefore conclude that non- 
linearities and non-local contributions have less impact on 
the third-order cross-statistics than on the corresponding 
auto-statistics. 

In order to measure the linear and the quadratic bias 
parameters we combine equations (ESI) and (I36II in a way 
which allows us to measure linear and quadratic contribu¬ 
tions (local or non-local) to the bias model independently 
from each other 


7 X _ 

br — 


S3 - Cl2 

qX 

•^3,h ~ ^12,1 

SsC^.h - Ci2S., 


r 

3^ 


(38) 

(39) 


This type of S 3 and C 12 combinations has been s uggested 
for th e au to-skewness and the aut o-correlators by ISzapudil 
lll998ll and Bejjnul MjaxmoiiH ll2012h and was further studied 
by Hoffmann et al.l'(2015bh. 


Sa.h 


S 3 + 3 ( c 2 



/bi 


(33) 


3.3.2 Bias measurements 


C\2,h — 


S. 


C 


C 12 + 2 (c 2 - - 32 ) 


/bi 


3,h — 


12,h — 


53 + C 2 - -g 2 


C 12 + C 2 — —g 2 


/bi 

/bi 


(34) 

(35) 

(36) 


(see Appendix lAj. On can notice the similarity with the 
corresponding equations for the three-point auto- and cross¬ 
correlations mu and ESI. The contribution of the second- 
order biases (local or non-local) to the cross-skewness and 
cross-correlators is three times smaller than in case of 
the auto-skewness and two time smaller than in the auto¬ 
correlator. 

In contrasts to the three-point correlations the second- 
order non-local contributions, described by equation El , can 
be taken into account as an effective second-order local bias, 

eff 2 t'o- 7 \ 

Cg = C 2 - -C/ 2 , (37) 

which consists of a local ( 02 ) and a non-local ( 32 ) contri¬ 
bution. This absorption of non-local contributions by 
results from the isotropy of spherically aver aged quant¬ 
ities, such as the skewness or the correlator (IChan et al.l 
I 2 OII already presented this effective local description of 
non-local biasing for spherically symmetric matter perturb¬ 
ations) . Hence, we do not expect the estimation of the linear 
bias via the equations (I33I) - (I36I) to be significantly affected 
by non-local contributions to the bias model. However, such 
non-local contributions will modify systematically the es¬ 
timation of the quadratic bias parameter C 2 . We will discuss 
the impact on the bias estimations in Section |4] 

By varying the bias parameters in equation (ISSll and 
EH) we can fit the measurements of the auto-skewness 
and auto-correlator for matter to the measured halo-matter 
cross-skewness and cross-correlator. We see in Fig. [5] that 
these fits are in reasonable agreement with the measure¬ 
ments for R > 20 /i“^Mpc. At smaller scales the hts lie 


Our measurements for the bias parameters, obtained from 
equations (1381) and (1391) . are displayed in Fig.[TO]with respect 
to the smoothing radius R. As for the measurements of C 12 
we set the distance between two spherical grid cells with ra¬ 
dius R to r = 2R. The measurements were performed for 
the four mass bins M0-M3 at redshift 0.0 and 0.5. We com¬ 
pare our results with corresponding measurements from the 
a uto-skewness and the aut o-correlator, and Cr, presented 
in iHoffmann et ^ ll2015bl l. The upper panel shows the lin¬ 
ear bias measurements together with the reference measure¬ 
ments from the two-point cross-correlation (equation (| 8 ))). 
which we consider to be a robust estimate of the true linear 
bias (see Section ED)- We hnd that, for all masses, 6(5 is 
less sensitive than br to the considered scale. This hnding 
is consistent with our previous conclusion from Fig. |9]that 
cross-statistics must be less affected by non linearities. The 
scale dependence is stronger for higher masses which can 
be attribu ted to non-linear and non-local terms in the bias 
functions (ISaito et al.ll20l3 l. This interpretation is suppor¬ 
ted by the fact that we find a better agreement between fej 
and fer at redshift 0.5 than at redshift 0.0. 

Our second-order bias estimations Ct and are shown 
in the bottom panel of Fig. HOI It is worth noticing that, 
as opposed to linear bias estimators, the estimator is 
affected by a larger scatter in the measurement than Cr. This 
can be explained by lower sensitivity of the cross-statistics 
{S^h ^r 2 h) non-linear and non-local contributions 
to the bias functions (i.e. cl®), compared the auto-statistics 
{S 3 ,h and Ci 2 ,h). This is indeed related to their monotonic 
behaviour with respect to the considered mass sample (see 
Fig. 0; the main dependence comes from the linear bias. 
The combination of and C ^2 h obtain therefore 
provides weak and noisy constrains on cl®. 

The fact that for highly biased tracers (M3) any non¬ 
local component in the bias relation is unable to explain 
why the two estimators br and b^ do not converge to the b^ 
estimator indicates the presence of a velocity bias, which 
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Figure 10. Top: Linear bias obtained from the b-r and estimators (respectively open and filled diamonds) compared to the reference 
linear bias from the two-point cross-correlation (dashed line) with respect to the smoothing scales R and at 2 : = 0.5 for the mass samples 
M0-M3 (Table Bottom: Quadratic reduced bias parameter C 2 estimated with Cr estimator (diamonds) at 2 : = 0.5. Error bars denote 
1 — cr uncertainties. 


has been shown to appear for very massive haloes (see 
IChan et al.ll20l'3) . This kind of contribution would add a 
dipole to the biasing relation which could modify signific¬ 
antly the br and estimators. 


4 RESULTS 

In the previous section we studied measurements of the lin¬ 
ear, quadratic and non-local bias parameters (for, C 2 and g 2 
respectively). These measurements were derived from three 
different methods which are all based on third-order statist¬ 
ics of halo- and matter density fluctuations. 

Two of these methods employ three-point auto-and 
cross-correlations. The first method is to compare the three- 
point cross-correlation with the three-point matter auto¬ 
correlation to derive the linear and quadratic bias paramet¬ 
ers via equation (HJ. This approach is based on the assump¬ 
tion of a local bias model. The linear and quadratic bias 
parameters, derived from this method are called fog and Cq 
respectively. The second method is to use particular combin¬ 
ations of three-point auto- and cross-correlations. The linear 
bias parameters, derived this way via equation (I23II . 6 aq, are 
independent of any quadratic contributions (local or non¬ 
local) to the bias function, as explained in the previous sec¬ 
tion. The quadratic and non-local bias parameters are ob¬ 
tained simultaneously by fitting predictions for the non-local 


component of the three-point correlation function, to 

::^Qcg, defined in equation (I24II . The quadratic parameters 
from such measurements is called caq- 

The third method is to use a combination of the halo- 
matter cross-correlator and cross-skewness with the corres¬ 
ponding auto-statistics for the matter field. The linear bias, 
derived from this so-called estimator via equation (I38II 
is called . The effective quadratic bias parameter, ob¬ 
tained from equation (I39II is called G ■ We demonstrated In 
the previous section that this effective quadratic bias (c^*^) 
consists of the quadratic bias and a non-local contribution, 
which cannot be distinguished from each other with the r 
method (see equation I37II . This feature is a consequence of 
the isotropy of the skewness and the correlator. Since we 
have shown that cross-statistics are weakly sensitive to the 
effective second order bias and therefore more sensitive to 
shot-noise we will use instead the Ct estimator. For clarity, 
in Table [2] we summarize our estimators and the notations 
that we use. 

In this section, we aim at comparing the bias estima¬ 
tions coming from these three methods. We therefore present 
the different linear and quadratic bias estimations for the 
four mass samples M0-M3 at the redshifts 0.0 and 0.5 versus 
the mean halo mass in each sample in Fig. 1111 Our bias meas¬ 
urements based on thee-point correlations are done using tri¬ 
angles with ri 2 = r' 13/2 = 36 /i“^Mpc. The estimations 























































Non-local bias contribution to third-order galaxy correlations 15 


Table 2. Bias estimators, definitions and notations 


Symbol 

Definition 

Equation 



Eq.® 


Qh — {Qm + Cq ]■ 

Eg.lllSIl 



Eq.iTT^ 

b^Q 


Eq.|[211 

b-r 

T 

Hoffmann et al. f2015b') 


'^h 


b? 

r ^ 

X 

Eq.|[381| 

CQ 

Qh = ^{Qm + cq} 

Eg.llTSll 


Qh=l^{Qrn+^-§-} 

Eq.dl 

CAQ 

CAQ = bi^lQh — Q^] 

Eq.l|25]) 

Ct 

S3^12,h—(^12S3^h 

'Th 

Hoffmann et al. ('2ni5b') 

C? 


Eg.IMt 


are based on fits of and between 26 7? ^ 40 /i“^Mpc 

using ri 2 = 2R configurations. 

The linear bias estimations from the different methods 
are presented in the upper panel of Fig. 1111 We compare 
these estimations to reference measurements from the two- 
point cross-correlation, defined in equation (1H]|, which we 
consider as reliable (see Section r3.ll) . The relative deviations 
to this reference linear bias are shown in the central panel. 

We find that the estimator bg, which neglects the non¬ 
local bias, overestimates the linear bias by 5-10%. The fact 
th at we found a s t ronger overestimation for 6 q (10-30%) 
in iHoffmann et al.l (l2015bl i can be attributed to the lower 
impact of non-local contributions to the three-point cross¬ 
correlation compared to the corresponding auto-correlation, 
as discussed in the previous section. 

The linear bias parameters from AQ is in excellent 
agreement with the reference for all mass ranges and at both 
redshifts. Deviations are in the range of the la of ^aq, while 
the latter roughly correspond to 1% of the amplitude. For 
the mass sample M3 deviations become slightly larger and 
more significant. We find this agreement also for smaller tri¬ 
angle scales, as we demonstrated for 2 = 0.5 in Fig. [7| and 
discussed in the previous section. 

Our linear bias measurements from the Tx estimator 
differs by less than 5% from the reference linear bias for the 
mass samples M0-M2 at both redshifts. Such deviation are 
comparable with the la errors and are therefore not signific¬ 
ant. At the highest mass sample M3 the linear bias derived 
from deviates signi ficantly from the referen ce. We found 
a similar behaviour in IHoffmann et ^ (l2015bli for the r es¬ 
timator. If this deviation was due to exclusion effect (leading 
to a wrong shot noise correction) then we would expect the 
estimator to be only weakly affected, which is not the 
case. An alternative explanation for these strong deviations 
can be the presence of velocity bias, which possibly cannot 
be neglected in this mass regime as we argued in Section r3.3l 

The lower panel of Fig. [11] shows how the methods 
compare in terms of estimating the second-order bias para¬ 
meter C 2 . Moreover, given that c-r only estimates the effect¬ 
ive second-order bias we also show the cl® (green triangles) 
computed from C 2 and g 2 (equation (I37ll i given by the AQ 
method. At redshift 2 = 0.0 we can see that the contribution 
of the non-local bias to the cl® is very small which, given the 


error on the c|® measured from Cr, does not allow to separ¬ 
ate the two measurements. At least, we see that for low mass 
bins the two effective second-order bias values agree within 
the errors. Hence, we can conclude that at redshift 2 = 0.0 
the effective second-order bias is a good approximation of 
the second-order bias. Regarding low mass bins at redshift 
2 = 0.5 one can see that, even if the effective C2 from AQ 
differs significantly from C 2 , it is difficult to say which one is 
in better agreement with c-r. The latter seems to be in better 
agreement with C 2 than the expected c|® coming from the 
three-point correlations. In case of the high mass bin, we 
observe again at both redshifts that the Ct estimator deliv¬ 
ers unreasonable results. We conclude that the Ct estimator 
can be used to estimate the second-order local bias para¬ 
meter as long as the considered tracer are not too massive 
(> 10^*h~^M q). Finally, we can say that the second-order 
bias, estimated from the three-point cross-correlation when 
the non-local bias is neglected leads to significant departure 
from the one obtained from AQcg- 


5 SUMMARY AND CONCLUSION 


We studied linear, quadratic and, for the first time in con¬ 
figuration space, non-local bias of halo clustering with re¬ 
spect to the clustering of the dark matter field. We there¬ 
fore employed various second- and third-order statistics of 
halo and matter density fluctuations in the MICE-GC sim¬ 
ulation. Our goal was to find if the overestimation of the 
linear bias parameter by the three-point auto-correla tion, 
whi ch we found previously i n Hoflhnann_^^_akU201^) (see 


also iManera and Gaztanagal I 2 OIII : IPollack et al.ll20121 , 
be attributed to shortcomings of the local quadratic bias 
model. Understanding this difference is crucial for breaking 
the degeneracy between growth and bias with three-point 
correlations, which would strongly amplify the statistical 
power of large-scale structure surveys. To achieve this goal 
we employ auto- and cross-statistics to disentangle the ef¬ 
fects of linear bias on second- and third-order halo statistics 
from those originating from non-linear and non-local bias. 

We started our analysis verifying how well the second- 
order clustering of haloes can be described by a linear bias 
model in Section 13.11 Comparing the amplitudes of the 
variances cr to those of two-point correlations ^ we found 
the former to be significantly affected by non-linear and 
possibly non-local contributions to the bias function. How¬ 
ever, the halo two-point correlations are well described by 
a linear local bias model down to scales of 10 /i~^Mpc. 
For this reason we employed such measurements from the 
two-point cross-correlation as an estimator for the linear 
bias, used as reference in the subsequent analysis. Note 
that we also expect departu res from scale independent bias 
(see IPesiacaues et aO l2010l 'l for highly bi ased samples at 
scales close the BAG peak (> 80 /i“^Mpc, IHoffmann et al.l 
l2015all^ . while the large errors at such scales should prevent 
a strong impact of this effect. 

For studying the impact of non-linear and non-local 
bias on the three-point correlation we compared in Section 
im bias measurements from the (reduced) three-point halo- 
matter cross-correlation to those from the auto-correlation 
from our previous study, using the local quadratic model. 
We found the linear bias from the cross-correlation to be 
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Figure 11. Top: Summary of the linear bias measurements presented in this work. Results are shown for the redshifts 2 = 0.0 and 
2 : = 0.5 (left and right respectively) versus the mean halo mass of each mass sample M0-M3, defined in Tabled (lower and upper limits 
of the mass samples are marked by vertical grey dashed lines. Symbols are slightly shifted along the the mass axis for clarity). The 
smaller panel shows the relative deviation of each estimator with respect to the linear bias from the two-point cross-correlation {h ^, 
solid lines). The estimators fog, b^q and (from equations 1I19II , 1I23II and II38I1 respectively) are displayed respectively by open circles, 
open triangles and open squares. For measuring 6g we assumed a local bias model, i.e. 52 = 0. Bottom: Summary of the second-order 
bias measurements Cg, caq and Cr from the same bias estimation methods, as used for the linear bias in the top panel (equation (|19|l h 
(fM and 1 I 39 II respectively). The Cj® (equation Il37ll i. obtained by combining C 2 and 32 from the AQ method, is represented with filled 
triangles. For estimations based on Q we used triangles with (ri2, ria) = (36, 72) /i“^Mpc configurations. Estimations using the t method 
are fits between 26 ^ iJ 40 /r“^Mpc. Error bars denote Icr uncertainties. 


closer to the reference than the linear bias from the auto¬ 
correlation. This is expected from second-order perturbation 
theory, which predicts the three-point cross-correlation to be 
less affected by quadratic local and non-local bias than the 
corresponding auto-correlation (compare equation 1)181) and 
GH)). However, the three-point cross-correlations delivers, 
as the corresponding auto-correlation, linear bias measure¬ 
ments, which lie significantly above the reference from the 
two-point correlation. 

To further verify if this overestimation can be attrib¬ 
uted to non-linear and non-local contributions to the bias 
model we take advantage of the fact that three-point auto- 
and cross-correlations are affected differently by non-linear 
and non-local bias, but equally by the linear bias (see again 
equation ca and IH). This property allows for combina¬ 
tions of the auto- and cross-statistics which isolate the linear 


from the non-linear and the non-local bias. We find the linear 
bias, measured by such a combination of three-point correla¬ 
tions independently of quadratic or non-local bias (equation 
1231) to be in excellent agreement with the reference from 
the two-point correlation (Fig. |T|. This finding is a strong 
indication that non-local terms are indeed the reason for 
the overestimation of linear bias from three-point correla¬ 
tion, when ignoring them by assuming a local quadratic bias 
model. This approach could be used to measure linear bias 
by cross-correlating galaxy with lensing maps. The presence 
of non-local bias also becomes apparent in our measurements 
of non-linear bias contributions (local and non-local) via 
equation (1241) . which are in good agreement with predictions 
for the non-local contributions to the three-point correlation 
(Fig. [ 5 ]). Our results therefore constitute the first detection 
of non-local bias in configuration space and demonstrate the 
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paramount importance of taking it into account when ana¬ 
lysing galaxy surveys. 


When the considered scales are two small (ri 2 < 
30 h“^Mpc), the non-local bias parameter 72, derived from 
our measurements, shows a strong scale-dependence, indic¬ 
ating the presence of higher-order local or non-local terms in 
the bias function. Instead, for scales larger than 36 /i“^Mpc 
we find a linear relation between the non-local and the lin¬ 
ear bias, over the whole mass range, as predicted by the 
local Lagrangian bias model (Fig. [8]). However, the amp¬ 
litude of this relation lies significantly below the local Lag- 
rangian prediction. This is in agreement w ith results from 
iBaldauf et al.l (l2012l f and lSaito_e^_alJ ((201^) , but in contra¬ 
diction with results from IChan et al] ( 2012l ~l. who find the 
non-local bias to be above the local Lagrangian prediction. 
Whether th i s latt er disagreement comes from the fact that 
IChan et al.l (l2012l ) analysed the bispectrum in Fourier space 
using different simulations is the subject of current investig¬ 
ations. An alternative reason for this discrepancy could be 
the inaccuracy of the prediction of the non-local contribu¬ 
tion to the three-point correlation of matter, from which we 
derive the non-local bias parameter. 


Besides three-point correlations we studied in Section 
bias from one- and two-point third-order statistics (i.e. 
the skewness and the reduced correlator respectively). These 
two statistics can be com bined into the so-called r-estimator 
dBel and Marinonill2012l f . which allows us to measure the lin¬ 
ear and quadratic bias parameters (the latter being local and 
non-local) independently from each other. An important dif¬ 
ference to the three-point correlation is that the r-estimator 
is an isotropic quantity. Non-local contributions to the bias 
function are therefore absorbed in an effective quadratic bias 
parameter (equation (1371) 1. Hence, the quadratic and non¬ 
local bias cannot be distinguished from each other by the 
r-estimator. Our measurements show that the r method de¬ 
livers a linear bias estimation which agrees at the la level 
with the reference linear bias from the two-point correlation 
at scales larger than 20/i“^Mpc, while the errors on the bias 
correspond to ~ 5% of the amplitude (Fig. 1101) . 

This result lines up with our findings from the com¬ 
bination of three-point auto- and cross correlations in the 
sense that third-order statistics is able to deliver accurate 
estimations of the linear bias when the employed estimator 
is independent of non-local and possible higher-order con¬ 
tributions to the bias function. 


Interestingly the linear as well as the quadratic bias es¬ 
timations derived from the r method are failing when meas¬ 
urements are performed in samples of very massive haloes 
{M > 10^“^ Mq). Since this effect appears in both, the auto- 
and the cross-statistics it might be attributed to velocity 
bi as, rather than non-Pois sonian shot-noise as we speculated 
in iHoffmann et ahl (l2015lJ l: further investigations are needed 
to understand this effect. 

Summarising our various bias estimations in Fig. [3 we 
find an overall variation in the linear bias of 10% with re¬ 
spect to the reference from the two-point cross-correlation. 
Deviations from the reference are more significant for meas¬ 
urements from the three-point correlation, which are based 
on the quadratic local bias model. The bias estimators, 
which are independent from non-local contributions {AQ 
and r) do not differ significantly from the reference for the 


linear bias, except for r in the high mass range, as we saw 
previously. 

Comparing the quadratic bias from the different es¬ 
timators we found the effective C2, obtained from the r- 
estimator, to be close to C 2 measured with the combination 
of three-point auto- and cross-correlations, which takes into 
account non-local bias. This finding shows that one can reas¬ 
onably neglect the non-local bias contribution when using 
the r-estimator to measure the quadratic bias parameter. 
This neglection is not possible for three-point correlations, 
as those deliver quadratic bias parameters which are signific¬ 
antly higher than measurements, which take non-local bias 
into account. 

Our results show that the local quadratic bias model 
is inadequate to describe halo bias in the MICE-GC simu¬ 
lation. Non-local second-order terms need to be taken into 
account for accurate measurements of the linear bias with 
three-point correlation function. Two approaches are pos¬ 
sible to do so. Non-local bias can be isolated from linear bias 
by combining different third-order statistics (i.e. r or AQ), 
or the non-local contributions need to be directly modelled. 
The first approach, on which we focused in this analysis, 
might be implemented in terms of cross-correlations between 
lensing and galaxy maps. We will test the second approach 
in a future analysis, but already provide here an expression 
for the non-local contribution to the three-point correlation 
in configuration space. At scales below 30/i“^Mpc we find 
indications for the presence of higher-order terms in the bias 
function (local or non-local). Modelling the non-local bias as 
a linear function of the linear bias parameter, as suggested 
by the local Lagrangian bias model, therefore appears to be 
only suitable at very large scales. 

In a future work (Bel, Hoffmann, Gaztanaga, in pre¬ 
paration) we plan to compare the linear and second-order 
bias measurement, obtained in the present analysis t o peak - 
background split predictions from I Hoffmann et al.1 (l2015ah 
and also bias parameters obtained by directly comparing 
halo- and matter over densities in the MIGE-GC simula¬ 
tion. These various bias estimations can be used to verify the 
universal relation between lin ear and non-linear bias para- 
meters, which we presented in iHoffmann et ^ ll2015al'l . 
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APPENDIX A: NON-LOCAL BIAS: A 
PERTURBATIVE APPROACH 

In this section, the Fourier transform /(q) or /(k) of a 
configuration space quantity /(x) or /(r) is defined as 
/(q) = (27r)-3//(x)e-‘"‘id3x. 

We assume that for large smoothing scales matter per¬ 
turbations are small such that they are well described by 
perturbation theory. We further assume that on such scales, 
the biasing function is described by the equation of this 
paper, which is 

:5.(x) = 6i |j^(x) + |(5^(x) - {5l)) + Je 2 (x)|, 

(Al) 

where 

52 ( x ) = - j Pi20v{(\i)Ov{^2)W\qi2R]e'^'^'^'^(Ecii(Eq2 

and /3 i 2 = 1 — ^ ^ i hence {Q 2 ) = 0. The velocity field 

V is normalised in such a way that at linear order its diver¬ 
gence is given by the Poisson equation = Sl- 

In the present calculation and in all this paper we work 
with density fields which are convolved with a Top-hat win¬ 
dow function of typical size R. In addition, we assume that 
perturbations are small and can therefore be described at 
second order by 

(5m(x) = (5l(x)-I- (5)^^(x), (A2) 

where Jn(x) is the linear contribution t o the fluctuations. 
And the second-order contribution (see iBernardeau et ahl 
I2OO2II can be expressed as 
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= f -F2(qi, q2)<5L(qi)(5L(q2)W[gi2-R]e*‘’^^ "‘d^qid®q2, 

(A3) 

where the second-order perturbation theory kernel is defined 
as F 2 (qi,q 2 ) = 1 -h (iT + S) “ 9 ^ 12 - This ker¬ 

nel shows that at second order, non-linearities arise due 
to mode-coupling. The factor g is the second-order growth 
factor which reduces to 2/7 in an Einstein-de Sitter uni¬ 
verse or as long as the growth rate of structures is given 
by / = . Th is function of time, c an be identified 

with the function B in ICatelan et al.l (Il995l ). or it is related 
to the function k (k + q = 1 / 2 ) introduced bv iBouchet et all 
lll992ll . a nd is aslo related to the u fun c tion (/j + 2g = 1) 
defined in iKamionkowski and Buchalteij (Il999l f. 

In the following, in order to use lighter notations we 
define v = 5h/bi. Using equation eu and expansion (IA2I) 
we can express v at second order as 

i2 = Sl + ^[SI-{SI)]+u‘-^\ (A4) 

where the quantity is arising from an effective mode 
coupling kernel such that 


Jl23 — 2(1 — 2g/3)^f2Cl3 + ^fl(’'i2,ri3) 

-|-i3ij(ri2, ri3)Li {^, 23 ) 

[C'H(ri2, ris) + 4:gDR{ri2, ris)] L2 (a‘23), (A8) 

where Ar, Br, and Cr and Dr are obtained from integ¬ 
rals over the linear power spectrum and gij = cosaij. In 
practice Bel (in preparation) has shown that in the large 
separation limit these functions can be identified with the 
fu nctions £(ri 2 ) (£12 in our notati ons) and (f){r 12 ) introduced 
bv iBarriga and Gaztana^ (120021') . It reads. 


Aij(ri2,ri3) 
Bfl(ri2,ri3) 
Ci?(ri2, ri3) 
DR{ri2,ri3) 


~ 0 

~ 0 
~ 1 
“ 3 


[^12 <!>'(rw) + ^i3(p'{ri2)^ 


+ 3 ^ 

ri2 


+ 3 ^ 

ri3 


(A9) 

(AlO) 

(All) 

A12) 


where </>(r) = J and ' = 

In addition, since we saw that g —>■ g + 'ii/b\ by substi¬ 
tuting it in equation (IA 8 I) one can split 77i23 in two contri¬ 
butions 


( 2 ( 2 )= f E|®(qi,q2)(5L(qi)(Iz,(q2)IT[qi2i?]e“’^2 "'d®qid®q2. 

(A5) 

Note that the effective second-order kernel is ob¬ 

tained by substituting g hy g + 72/&1 in F 2 . From equa¬ 
tion (pH) . taken at three positions ri, r 2 and rs, we can 
express the three-point correlation function of u, = 
(i/(xi)i/(x 2 )i^(x 3 )) we obtain 


C = + ^[{Sl{ri)SL{r2)SL{r3))+perm] 

[^23(<5i) -I-perm] 

+(! 2 ^^'(ri)( 5 i(r 2 )( 5 L(r 3 )) + perm, (A 6 ) 

where use the definitions = {5L{ri)SL{rj)) and = 
(5i:,(xi)5r (x 2 )dn(x 3 )) . Assuming that the linear part of the 
density field (5 _l is Gaussian, then its three-point correlation 
function is n ull and its f our-point expectation value can be 
expressed as (lFrvlll98'4^ 

((5r(ri)5i,(r2)(5r(r3)(5i,(r4)) = ^12^34 + Ci3^24 + ^14?23 

from which we can derive (ri = r 4 ) that 

(di(ri) 5 i(r 2 )< 5 L(r 3 )> = + 2^Ut3- The first 

term cancels with the third term of equation (IA 6 I) . As a 
result we obtain that 


= C 2 |]Ci 2 ^i 3 -I- perm] -I- J123 + perm, (A7) 

where J 123 = (i 2 ^^^(ri) 5 _L(r 2 )( 5 z,(r 3 )). This expression has 
been calculated by Bel (in preparation) who has shown that 
taking into account the convolution with the spherical Top- 
hat filter is very well approximated by a Legendre expan¬ 
sion at second order (which turns out to be exact when no 
smoothing is applied). Hence 


v7i23 = Fi 23 -|-p2Afl23, (^13) 

where /C 123 = | [ri23 —^ 12 ^ 13 ] is the non-local part which 
has been introduced in equation dsni) and ri 23 can be ex¬ 
pressed as 


fl23 = 2(1 - 25r/3)^f2^13 + ^fl(’'12, 7’13) 

+BR(ri2,ri3)Li(g,23) 

+ [C'fl(ri2, ri3) -I- 4/3gDR(ri2,ri3)] L2(m23)(A14) 

From equation (IA14I) , we can express the three-point correl¬ 
ation function of the matter field as 


Q 


FT 

m 


1 -I- r23i -i- r3i2 


C 


m 

H 


(A15) 


which is equivalent to the expression given in 

iBarriga and Gaztanagal (120021 '). By analogy, we can 
define a non-local three-point function Qnioc as 


Qnloc — 


/C123 + /C23I + /C312 

/-m 

Sjf 


(AI 6 ) 


The definition (IAI 6 II can be written in terms of the function 
ri 23 and its permutations 


^ 2 f Fi 23-l-r23i-l-r3i2 ) ^ ^ 

Qnloa=^<^ - ^ -Ij, (A17) 

where ri 23 = DR{r\ 2 ,ri 3 )L 2 {g 23 ) which contains only a 
quadrupole contribution. On the other hand, from equation 
(IA 6 I) . we can express the reduced three-point function of 
haloes Qh = 6 r^C'"/Cir with respect to ri 23 and /C 123 




(A18) 
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Finally, using definitions introduced before, we find that the 
reduced three-point function of haloes can be expressed as 


Qh — {Qm -|- C2 -|- g2Qnloc} ■ (^^ 19 ) 


In the end, from equation (IA19I) we can take the limiting 
cases which correspond to the skewness and the correlator. 
In case of the skewness Bel (in preparation) has shown that 
only the monopole contributes. As a result ri 23 and all its 
permutations are null therefore 



(A20) 


The second limiting case is obtained when ri 2 = ris and 
ct23 = 0 , in this case the permutations r23i and r3i2 are 
null so only the ri 23 is contributing. As a result we get 



(A21) 


however the function I?jj(ri 2 , ri 2 ) contributes by less than 
5%/2% on scales greater/below than R — 20h~^Mpc com¬ 
pared to Ci 2 . As a result considering that this term is mul¬ 
tiplied by 2/3g2 it will contribute in total at most to 2% 
(high mass bin) and to less than 0.1% for lower mass bins 
and can be therefore safely neglected, it reads 



(A22) 


The same kind of approach can be used to generalize the 
relations IA191 IA20I and IA22I to their corresponding cross¬ 
statistics. Note that those relations are valid at the three- 
level and do not imply the large separation limit approxim¬ 
ation. 










